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We study the quantum Goos-Hanchen(GH) effect for wave-packet dynamics at a nor- 
mal/superconductor (NS) interface. We find that the effect is amplified by a factor (Ef/A), with 
Ef the Fermi energy and A the gap. Interestingly, the GH effect appears only as a time delay 
6t without any lateral shift, and the corresponding delay length is about {Ef / A)Xf , with \f the 
Fermi wavelength. This makes the NS interface "sticky" when A <C -Ef, since typically GH effects 
are of wavelength order. This "sticky" behavior can be further enhanced by a resonance mode in 
NSNS interface. Finally, for a large A, the resonance-mode effect makes a transition from Andreev 
to the specular electron reflection as the width of the sandwiched superconductor is reduced. 

PACS numbers: 03.65.Vf, 74.45.-|-c, 73.40.-c 



Interference is an important phenomenon in both 
quantum and wave mechanics, and is often responsible 
for intriguing effects not understandable from a classi- 
cal mechanics and ray dynamics. The Goos-Hanchen 
(GH) effect [l|, first predicted by Newton, is one ex- 
ample: When an optical beam is totally reflected at a 
dielectric interface, a lateral shift along the interface is 
induced. The GH shift has been studied in various quan- 
tum and wave systems, amongst which graphene 0, 
optical beams [j] , dielectric microcavities [5|-|7[ and neu- 
trons 

The GH shift in optics can be understood by the depen- 
dence of the phase loss occuring at total internal reflec- 
tion on the angle of incidence. An analogous shift is 
observed in the quantum mechanical step potential prob- 
lem: In the one dimensional case, if the energy E of an 
incident particle is less than the potential height V, the 
particle is totally reflected with reflection coefficient is 
R = e^'^ (with real tp). The phase (/? increases monoton- 
ically from — tt to zero as the energy E varies from zero 
to V . If we consider an incident Gaussian wave packet 
with velocity = Hko/m, 



with the time delay given by 



'ilji{x,t) = / dkg{k)e 



^ikx —iEt/h 



(1) 



where g{k) — e~^^~^^'^ , cr is the width of the Gaussian 
window, and E = (hk)^ /2m, the wave packet reflected off 
the step potential at a; = becomes 



(2) 



The stationary phase approximation gives a time-delayed 
motion of the centre of the wave packet, x = —vo{t — St), 



^' = ^dE- 



(3) 



As a result, the delay length lo = voSt is Z^i ~ Aq with 
Ao = 2-K/kQ the central wavelength of the wave packet. 
For a wave packet propagating in two dimensions with 
y-direction velocity voy, the time delay results in a lat- 
eral shift of Igh = voySt, analogously to the optical GH 
shift. Note that the GH shift here is of the order of a 
wavelength. This 6t is similar to the Wigner delay time 
in scattering problems [lo| . The phase change ip is also 
related to the non-integer Maslov index for the quantiza- 
tion condition [Tlj. 

On the other hand, the behavior of an electron or a 
hole at normal/superconductor (NS) interface has been 
studied extensively. For a small superconducting pairing 
gap A (^ lQ~'^Ep) an incident electron from the normal 
side is retro-reflected as a hole when the excitation en- 
ergy is less than the pair gap A of the superconductor; 
this is called Andreev reflection [l^l , and has been stud- 
ied in systems as diverse as graphene [l3j and bosonic 
condensates IJ] . It has been directly experimentally ob- 
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served [15|. In addition, the advance of laser cooling 
techniques for atoms allows experimental realization of a 
similar interface with a large A comparable to the Fermi 
energy Ep, where specular reflection of particles (in ad- 
dition to Andreev reflection of holes) occurs with appre- 
ciable amplitude pli-[l8|. In spite of much attention to 
the electronic properties at the NS interface, the phase 
change under reflection has not been explored so far. It 
is therefore desirable to understand the influence of the 
phase change on reflection off the NS interface and the 
associated GH effect. 

In this Letter, we study reflection off NS and NSNS 
interfaces and find that the GH effect is amplified by a 
factor {Ep/A) compared to that of other interfaces. This 
is because at the NS interface the incident electron can 
be reflected within the energy scale of A, while the wave- 
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FIG. 1: (color online) (a) Schematic diagram for electron and 
hole excitations near the Fermi surface in (fc^,, ky) plane. The 
arrows indicate the direction of group velocity of the corre- 
sponding wave packets, (b) The real-space illustration of the 
excitations in (a). An incident electron can be bounced into 
a hole (A) and an electron (B). 1^ and ly^ are GH shifts ex- 
pected from a naive extension of the step-potential case, but 
they turn out to be zero in this NS interface (see Fig. [5] and 



length is determined by the whole energy scale of ^ Ep. 
Moreover, the GH effect appears as a time delay with- 
out the lateral shift which has been the hallmark of the 
GH effect in other interfaces. The corresponding delay 
length is then given by Id ^ [Ep / ^)^f , where Xp is the 
Fermi wavelength, implying that the NS interface is very 
"sticky" in comparison to typical wavelength-order delay 
lengths in optical and condensed-matter interfaces. We 
also show that the time delay can be increased by a reso- 
nance mode in NSNS interface. Finally, for a large gap A 
(~ Ep) the resonant reflection exhibits a transition from 
Andreev retroreflection to specular electron reflection as 
reducing the width of sandwiched superconductor. 

In what follows, we describe our theoretical setup, then 
begin our study with the A/Ep ^ 1 for both NS and 
NSNS interfaces. We then turn to the general A/Ep 
case, for which we discuss the crossover between Andreev 
and specular reflection before drawing our conclusions. 

In a superconductor, electrons and holes are coupled 
to each other with a coupling strength given by the pair 
gap A; this situation is described using the Bogoliubov- 
de Gennes (BdG) equations for a two-component wave 
function [ll,[2^: 



H A 

A -H 



(4) 



H 



2m 



(5) 



where e is the energy measured from the chemical po- 
tential fi (ii ~ Ep in our case) while u and v arc the 
electron and hole wave functions, respectively. The pair 
gap A is zero in a normal conductor, i.e., the electron 
and hole are independent of each other. We also assume 
that A is a real number, that is, we consider an s-wave 
superconductor and a static situation. 

Wc begin by focussing on the case of an incident elec- 
tron with an energy e ( < A) as schematically illustrated 
in Fig. [1] From the BdG equation, the two independent 
plane- wave solutions [ll,[2l[ with the same y-dependence 



are 



1 
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where — ^Je^ ~ A^ and 



2m , 
1? 



2m 
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where kpx = 
wavenumber kp when k 



2m 77' U2 

-j^Ep - ky 



and it becomes Fermi 

y — 0. Note that in the nor- 
mal conductor with zero A, the solutions become inde- 
pendent electron (u) and hole (v) wavefunctions since 
[e — Vl)/A = in Eq.(|ni), and their wavenumbers, k^ 
and k^ , are given by Eq.([7]) with D, = e. 

We obtain the reflection coefficients for the reflected 
electron and hole plane waves by matching wavefunctions 
at the interface (x = 0). The wavefunction in the normal 
part is 



N 



'ik' X 
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(8) 



while the wave function in the superconductor is ^5 with 
decaying u and v for our case of e < A. From ^'Ar(O) = 
^'5(0) and ^''^(0) = \E's(0), we can obtain the reflection 
coefficients A (hole) and B (electron). 

In order to study a dynamical situation, we now switch 
to an incident wave packet rather than a plane wave. 
The incident electron with an average wavevector fco ~ 
(kox, koy) can be described as a Gaussian wave packet. 



ui{x,t) — / dkxdkyQ{k)e 



i{k^x+kyy—et/h) 



(9) 



where g{k) = e^l'^-'^olV'^" and e = {h^\k\'^ /2m) - Ep. 
Note that at i = the wave packet hits the interface 
located at x = 0. The reflected hole is then represented 
by 



VFi{x,t) = I dkxdkyQ{k)Ae 



i{k^ x+kyy—Et/h) 



(10) 
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FIG. 2: (color online) Wave-packet dynamics when A — 
W^'^Ef- (a) The phase (pA (red line) and the Gaussian win- 
dow function g{e) (blue line) used. In the lower pannel, elec- 
tron wave packets (red line) at t = —8 and reflected hole 
wave packets (blue lines) at t = 6, 8 are shown, (b) The 
same as (a) for the resonant reflection in the NSNS inter- 
face with {d,L) = (1,4). (c,d,e) Wave-packet dynamics with 
45 degree incident angle. The scale of the normal conductor 
shown above is (100 x 200) fim^. (c) The incident-electron 
wave packet ai t — —10. The yellow arrow denotes the inci- 
dent path. The white arrows show schematically a typical GH 
shift of the step-potential case, (d) Andreev reflection. The 
reflected-hole wave packet at f = 10 is shown, (e) The same 
as (d) in the NSNS interface with (d, L) = (0.8,2.8). The 
delay length corresponds to the distance between initial (red 
dot) and final position (white dot). Here, the units of length 
and t is /xm and 6.58 x 10~^^s, respectively, (see supplemental 
material for wave-packet animations) 

and the wave packet for the reflected electron can be ob- 
tained by the replacements A ^ B and fc^ — >■ — fc^: in the 
above equation (see Fig. [T](a)). For the normal incidence 
case koy = 0, it is easy to get time delay from the sta- 
tionary phase point: The peak position of the reflected 
hole is given by 

where ipA is the phase of A. Similarly, for the reflected 
electron wc find 5tf^ = h{d(pB/de). 

We now concentrate on the case of relevance to solid- 
state and liquid helium experiments [iHl, A ^ Ep, 
specifically with A/Ep = 10~^ which is the typical value 
in solid-state superconductors. We set the Fermi energy 
as Ep = 5.5 eV throughout this Letter. This case corre- 
sponds to Andreev reflection [l^ , so that only the retro- 
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FIG. 3: (color online) (a) \A\'^ (red) and \B\^ (blue) are 
shown as a function of e, when A — 0.2Ef (dashed line) 
and A = 0.8i5F (solid line), (b) The corresponding phases, 
ifiA (red) and (ps (blue), (c) Change of \Af (red) and 
(blue) evaluated at e = 0.5A with increasing A. (d) The 
corresponding delay time 5t^ as a function of A. 

reflected hole has appreciable amplitude (|^| ~ 1 and 
\B\ ~ 0). The phase cpA is drawn in Fig. ^ (a) and it 
is well approximated by ipA = arccos(e/A) when B = 0. 
This phase variation enables us to calculate the delay 
time and length, defined by Z^i = vpSt-^, vp is the Fermi 
velocity. We emphasize that the delay length is amplified 
by a factor (Ep/A), i.e.. 

Id ^ (5^) A^, (12) 

which is a distinctive feature of the NS interface, not- 
ing that typical shifts in other interfaces arc wavelength- 
order. Thus the NS interface can be characterized as 
"sticky" due to the long-time stay at the interface. 

Since the delay time is proportional to the slope of 
the phase change, dipA{B)/d£, it can be extended by in- 
creasing the slope. This may be done with the help of a 
resonance mode introduced by changing the NS interface 
into a NSNS structure. We denote the widths of first 
superconductor barrier and next normal conductor as d 
and L, respectively. The phase (fAi^) for plane waves in 
the NSNS case is shown in Fig. [2](b): the slope becomes 
steeper as d increases, eventually approaching an abrupt 
2tt phase jump at the resonance point. 

Now consider a Gaussian wavepacket, beginning with 
the normally incident case, ko = {kox, 0). For the window 
fimction g{k) in Eq. (O, we take /cqx = , where 

kM = ^jkl + |?A, and a = , n = 4 and 5 for the 

NS and NSNS interface. The k^ dependence of the win- 
dow function is given by g{kx) = g{kx,0), and it corre- 
sponds to g{e) = e"(''"^o'^/'^s Eq = A/2 and cr^ = A/n, 
which are shown as blue line in the upper panncls in [2] 
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FIG. 4: (color online) Various resonant reflections when A = 
0.2Ef. (a) [d, L) = (3, 13). (b) (d, L) = (4.7, 13). (c) (d, L) = 
(20, 13) in the unit of A. The red dots denote the initial 
position and white and yellow dots correspond to the final 
positions of hole and electron, respectively. Note that the 
incident-electron velocity is not parallel to the reflected-hole 
velocity in this large A case, because both are parallel to 
radial lines of Fermi circle as shown in Fig. [T] (a). The scale 
of the normal conductor shown above is (0.12/im x 0.24/im). 
(d),(f),(g) \A\'^ and \Bf as a function of e/A for the cases 
of (a),(b),(c), respectively, (e) The phases, (pA and ipB for 
the (d, L) = (3, 13) case (see supplemental material for wave- 
packet animations). 



(a),(b). The wave-packet movements are shown in the 
lower pannels. The red line and blue line describe the 
incident-electron and reflected-hole wave packets, respec- 
tively. For the NS interface in Fig. [2] (a) , we can see some 
distance between wave packets of electron at t = — 8 and 
hole at i = 8, correspondig to the delay length, In = vpSt 
(note \ve\ — \vh\ ~ vp in this case, Ve and Vh are electron 
and hole group velocities). For the resonant reflection 
in Fig. [2] (b) , a longer time delay is expected from the 
steeper slope of fAie), and the extended delay length is 
clearly shown in the lower panel of Fig. [2] (b). 

Now consider a finite incidence angle, i.e., koy ^ 0; 
here one might expect a GH shift for the reflected hole 
like the Ih in Fig. [T] (b), based on the separablility of 
X and y coordinates like for a diagonal step potential. 
However, this is not the case. The wave-packet dynamics 
are illustrated in Fig. [2](c)-(e). Figure [2] (c) shows the 
incident wave packet with 45 degree incident angle. The 
retro-reflected hole packets are shown in Fig. [2] for the 
NS interface (d) and for the NSNS interface (e), for two 



times symmetrically before and after the time of impact 
on the interface. From the distance Id between initial 
and flnal position of the wave packet it is clear that there 
is a time delay, but no GH shift at all. This can be un- 
derstood by the dependence of the energy e on ky, i.e., 

^ ~ Wni^x + ^y) ^ ^F- This glvcs risc to a time delay 
Sty in the y-direction too, and this time delay is the same 
as 5t^ . Thus we see only time delay without GH shift 
in Fig. [2] (d),(e). The deformed shape of the hole wave 
packet comes from the relation between kx and k^ . In 
other words, the symmetric Gaussian window Q{kx,ky) 
becomes (approximately) an asymmetric Gaussian win- 
dows in {kj^,ky) plane, as depicted in Fig. [T](a). 

We now turn to a discussion of the more general case 
where A/Ep is not small. As A/Ep increases the nature 
of the reflected object changes from purely Andreev- to 
mixed Andreev- and specularly-reflected [16j . Figure [3] 
shows how the coefficients A, B and the delay time 6t 
change as A increases. \A\'^ and \B\'^ and the phases (pA 
and ifB are shown, as a function of e, in Fig. [3](a),(b), 
where the dashed lines are for A = 0.2Ep case and the 
sohd lines are for A = O.SEp case. When A = 0.2Ep, the 
cocfficents indicate almost pure Andreev reflection (|A| = 
1 and \B\ = 0). At A = O.SEp, \B\'^ becomes appreciable 
so specular reflection becomes important. Fig. [3] (b) 
implies that {d(pA,B /de) ^ 1/A so that the delay time 
behaves as 5t ^ 1/A, which is drawn in Fig. [3] (d). In 
Fig. [3] (c), I^P, l-Bp measured at e = 0.5A are shown 
as a function of A, indicating again the non- negligible 
specular electron reflection for a large A. 

It is interesting that the resonant reflection in the large 
A case exhibits a dramatic change of the reflected object 
depending on the width of the superconductor barrier, d. 
As an example, the various reflections, when A = 0.2Ep, 
are shown in Fig. U) Here we take L = ISA, which sup- 
ports a resonance mode around e = 0.5A for the 45 de- 
gree incident angle. As shown in Fig. [4] (a), both electron 
and hole are reflected out at d = 3 A, and then electron- 
dominant reflection is observed at d = 4.7A. contrast to 
Andreev reflection of the small A case. For d ~ 20A, we 
see the hole-dominant reflection, back to Andreev reflec- 
tion. 

This rather peculiar behavior can be understood as 
an effect of the resonance mode. Figures 5] (d-g) show 
variation of |Ap and around the resonant point for 
the above three different cases. For thin barrier case, 
d = 3A, the resonance mode is quite leaky so that the 
resonance mode affects on the coefficients within a rather 
broad range in e, resulting comparable values of \A\'^ and 
at the resonant point (see Fig. 111(d)). Their delay 
time has a similar value as expected from the slope of the 
phase variation in Fig. 2] (e). As increasing the barrier 
width, the broad range is getting narrower and reach its 
maximum, \A\'^ = and \B\^ = 1 (see Fig. g] (f)). If the 
barrier width is getting thicker, the single peak struc- 
ture starts to split into a double-peak structure, show- 
ing \A\'^ = and \B\'^ = 1 at each peak positions, and 
the width of peaks is getting narrower as shown in Fig. 
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|4] (g). This double-peak structure of a resonance mode 
is a charactericstic feature of two-component resonance 
where each component has its own wavelength [22j . 

In conclusion we have studied wave-packet dynamics at 
NS interface, and found that the GH effect is amplified by 
the factor Sf/A. Interestingly, the effect appears only 
as a large time delay, hence "sticky", without any lateral 
shift. We also demonstrated that in NSNS interface the 
GH effect is even further enhanced by a resonance mode 
and for a large A the transition from Andreev to specu- 
lar electron reflection occurs when the decay rate of the 



resonance increases. Although we have focussed on the 
incident electron case, all intriguing features discussed in 
this Letter can also be found in the incident hole case. 
We believe that the time-delay effect can be confirmed by 
an experimental observation, since direct observation of 
Andreev-reflected electrons in clean systems (where the 
motion is practically ballistic) is possible [T5| . 
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